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ABSTRACT 

Properties of low-energy (infra-red) quadrupole modes (LEQM) of multipolarity Xfi =20, 
21 and 22 in deformed sodium clusters are studied within the Kohn-Sham LDA RPA ap- 
proach. Possible manifestations of LEQM in stimulated Raman adiabatic passage (STIRAP) 
reaction are discussed. It is shown that, in free light clusters, where the low-energy spectrum 
is delute, LEQM can be unambiguously identified as particular electron-hole pairs. This gives 
a chance to reconstruct the mean field level scheme near the Fermi surface. Moreover, due 
to the connection with electric A// =21 mode, the scissors mode can be detected. In heavy 
(supported) oblate clusters, LEQM are in general rather involved. Nevertheless, some inter- 
esting A// =21 and 22 structures determined by specific deformation effects can be resolved. 
The origin of the structures is discussed in detail. 
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I. INTRODUCTION 



Collective oscillations of valence electrons in metal clusters manifest themselves in a 
variety of electric and magnetic (orbital) plasmons (see |l| for review). Till now only the 
electric dipole (El) plasmon was thoroughly investigated 0]. Other plasmons were not 
still observed and our knowledge about them is scares. At the same time, they promise 
interesting physics and so devote careful study. 

In the present paper, we consider low-energy electric quadrupole modes in axially de- 
formed clusters. Quadrupole excitations of valence electrons in such clusters are collected 
into bunches AN = and 2, where N is the principle shell quantum number P, 0]. The 
low-energy bunch AN = corresponds to E2 transitions inside the valence shell. It is clear 
that this bunch exists only in clusters with partly occupied valence shell. Just these clusters 
exhibit a non-spherical shape. The high-energy bunch AN = 2 corresponds to E2 transi- 
tion through two shells. This bunch exists in clusters of any shape and exhausts most of 
the quadrupole strength. Though the bunch AN = 2 shares the energy region with the 
dominant El plasmon, it can be possibly observed by means of angular-resolved electron 
energy-loss spectroscopy (AR-EELS) at electron scattering angles ~ 6° 

In this paper, we will consider the low-energy quadrupole modes (LEQM) with AN = 0, 
whose properties, by our knowledge, were never earlier investigated. We will show that 
LEQM demonstrate spectacular structures connected with cluster deformation. 

LEQM are probably too weak to be observed in AR-EELS but have a chance to be 
detected in two-photon processes like Raman scattering (RS) or resonance fluorescence (RF) 
where the competition with El modes is completely excluded. These processes start from 
excitation of an intermediate El electronic state (virtual or real). Then the coupling of the 
intermediate and low-energy states results in the population of the latter, again through 
El transitions (e.g., the electron-phonon coupling results in RS population of ionic phonon 
states) . 

As far as we know, there was only one attempt to observe LEQM by RS in nanoparticles 
. In that experiment a broad bump was observed at ~ 0.12 eV in supported silver clusters. 
The electronic nature of the bump was justified by the inverse cluster size dependence of the 
bump energy. Recent data of this group hint at the deformation dependence of the bump 
energy 0. Explanation of these data requests a strong theoretical support. Our study could 
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be a first guide in this new field. 

Perhaps, a weak coupling of electron states does not favor their detection by RS. In this 
connection, the stimulated Raman adiabatic passage (STIRAP) seems to be more appro- 
priate. This method uses partially overlapping pulses (from pump and Stokes lasers) to 
produce complete population transfer between two quantum states |1 > and |3 > through 
the intermediate state |2 > jl^. The combination of the pump and Stokes frequencies has 
to be resonant with the two-photon Raman transition. In fact, the Stokes pulse plays the 
same role as the coupling between the states |2 > and |3 > in RS. 

LEQM seem to be most interesting in the cases of free light non-spherical clusters and 
supported heavy oblate clusters. In free light clusters, the low-energy spectrum is very dilute 
and thus LEQM can be safely resolved. We will show that LEQM can be unambiguously 
identified as particular electron-hole (leh) pairs involving the single-particle electron states 
near the Fermi level. This allows to obtain, at least partly, the single-particle spectrum in 
the clusters. Since this spectrum is sensitive to many factors (deformation, temperature, 
influence of detailed ionic structure, ...), we thus handle the effective method to study diverse 
cluster's properties. The spectrum can be also used as a robust test for theoretical models. 



the latter 



Besides, due to the connection between LEQM and scissors mode 
can be detected as well. 

Heavy clusters in experiments are usually supported and range shapes from spherical to 
strongly oblate Size and shape of supported clusters can be monitored If heavy 
clusters are supported on dielectric (and porous) surfaces, then the interface effect is mainly 
reduced to shaping the cluster [2] . Such clusters can be simulated, in a first approximation, 
by free clusters at a certain oblate shape. In the present paper, we pretend only to a rough 
estimation of LEQM properties and so this simulation is reasonable for our aims. 

As a simplest case, we will consider sodium clusters. However, LEQM properties de- 
scribed below are of a general character and so should also take place in other clusters 
exhibiting quantum shells, e.g. in clusters of noble metals. 

II. CALCULATION SCHEME 

LEQM are described in the linear regime within the random-phase-approximation (RPA) 
method based on the Kohn-Sham LDA functional llJl . The ions are treated in the soft 
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jellium approximation. High reliability of the method has been checked in diverse studies of 
El plasmon in spherical ^3 ^^^1 deformed jj, Q| clusters. 

We consider axial prolate Na^/g as an example of a free light cluster. Equilibrium 
quadrupole and hexadecapole deformations of Naj'g (^2 =0.57 and =-0.18) are deter- 
mined by minimization of its total energy P, Q|. The neutral Nang at oblate deformations 
S2 =-0.3 and -0.5 is used as a simulation of heavy supported clusters. 

STIRAP is treated as a resonance fluorescence (RF) process running through A// = 10 
and 11 branches of the dipole plasmon, as intermediate states (see Fig. 1). In strictly 
deformed nuclei, where the branches are well separated by energy (see Figs. 2 and 3), 
independent measurements through every branch are preferable. The population picture is 
then simpler and allows easier interpretation. For example, if STIRAP runs only through 
A// = 10 plasmon branch, then the population of the low-energy quadrupole 22 mode is 
forbidden and thus remaining 20 and 21 modes can be easier treated. 

STIRAP population of LEQM is approximately calculated as a coherent sum of indepen- 
dent two-step processes, every one being a product of dipole photoabsorption and emission: 

^2M2i2 = J]'^^'it(0 ^ l/^l^l) ■ ^lu(l^l^l ^ 2^2^2) (1) 

h 

=^ X^^iMiiil < Ql^J.lil\El\0 > p(u;i^^ij 
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Here \lfiiii > and |2/i2«2 > are RPA states of the dipole plasmon and LEQM, respectively. 
Index i runs all the RPA states of the given multipolarity in the chosen energy interval. The 
operator of the dipole transition is electron-hole in the photoabsorption and electron-electron 

(hole-hole) in the photoemission. 

Eq. follows from the general expression for RS cross-section [17] if one neglects the 
interference between the neighbouring RPA states. This is the case for light clusters (see Fig. 
2) but not for heavy ones whose spectrum is rather dense jj, . However, the interference 
mainly leads to smoothing the response, which can be taken into account by the reasonable 
averaging the results. In the present paper, we weight the low-energy responses by the 
Lorentz function with the averaging parameter A =0.1 eV. This should simulate both the 
interference and temperature smoothing. We emphasize that the structures discussed below 
are strong enough to be completely smoothed out and so the main conclusions of the present 
study have to be valid in spite of the interference. 
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Spread of collective oscillations is known to increase with the excitation energy. So, for 
the dipole plasmon which lies much higher than LEQM, we use the larger averaging with 
A =0.25 eV. Such averaging was successfully used in our previous RPA calculations for the 
photoabsorption [4, 15, Ig j. 



III. RESULTS AND DISCUSSION 



Main results of the calculation are collected in Fig. 4. The first line of the figure contains 
photoabsorption for low-energy quadrupole states A^u =20, 21 and 22. Photoabsorption 
is a useful step in any analysis of electron modes. The second line of Fig. 4 exhibits 
photoabsorption for the scissors magnetic dipole mode. The next two lines provide the 
STIRAP populations for the cases when the reaction runs separately through A// = 10 
and AyU = 11 branches of the dipole plasmon (these cases are marked in Fig. 4 as RF ElO 
and RF Ell, respectively). 

In the energy intervals given in Fig. 4 (0-1.5 eV for Na^g and 0-1.0 eV for Nans), we take 
into account all the quadrupole RPA states. As for the dipole intermediate RPA states used 
in the STIRAP calculations, we involved all the states in Na^g and 40-50 of most collective 
states in Nans, determined by the appropriate cut-off of the photoabsorption strength. 



A. Light free clusters 

The first plot in the left column of Fig. 3 shows that LEQM photoabsorption in Na^'s is 
dominated by two A// =21 peaks. Following the level scheme for Na^g given in Fig. 5, these 
peaks can be associated with electron-hole pairs [211]- [202] and [211]- [200]. Our analysis 
shows that i) collective shifts for these peaks are very small and ii) just the above leh 
pairs provide dominant (up to 95%) contributions to the normalization conditions of the 
peaks. This justifies identification of the peaks as [211]-[202] and [211]-[200] leh pairs. The 
contributions into the optical response of other quadrupole modes, 20 and 21, given in Fig. 
5 is negligible because of low values of their transition single-particle matrix elements. 

Both [211]-[202] and [211]-[200] leh pairs include the Fermi level [211] whose energy can 
be obtained from the ionization potential data. Then one immediately gets the energies 
of electron levels [202] and [200]. So, STIRAP measurements together with other data 
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(ionization potential, photoemission data, etc) allow to obtain the single-particle spectrum 
in light deformed clusters. Since this spectrum is sensitive to many factors (deformation, 
temperature, influence of detailed ionic structure, ...), we can use it to study diverse cluster's 
properties. Beside, the spectrum can serve as a robust test for theoretical models. 

The next plot for Naf^ exhibits the photoabsorption cross section for the scissors mode. 
We see again two peaks corresponding to [211]- [202] and [2 11]- [200] states discussed above. 
While the lower state [211]-[202] favors the Ml response, the higher state [211]-[200] responds 
mainly to E21 field. Both photoabsorption plots display a close connection (coupling) be- 
tween the scissors Ml and quadrupole E21 modes. The modes are characterized by one and 
the same set of leh pairs with quantum numbers A'^ = 1+ (where A is the projection of the 
orbital moment to the symmetry axis z and n is the space parity of the state). The coupling 
of electric and magnetic modes with the same quantum numbers A'^ is a general feature of 
deformed quantum systems. For example, this feature is well known in atomic nuclei jigj ). 

In fact, both scissors Ml and quadrupole E21 modes represent one and the same intrinsic 
electron motion in a deformed cluster. The scissors mode is more specific, while the E21 
mode is more general. To illustrate this point, we expand the single-particle wave functions 
in a deformed mean field in terms of the spherical basis {nLA) 

^.=[Arn.A] = J2<LRnL{r)YLAin). (2) 

nL 

This allows to evaluate the single-particle orbital Ml transition amplitude between hole 
(z/ = h) and particle (z/ = p) states 

nL 

Eq. (j2)) shows that the scissors mode is generated by Ap = A/^ ± 1 transitions between 
the components of one and the same spherical (riL)-level. Such selectivity is explained by 
the fact that the scissors operators, and Ly, do not depend on the space coordinate 
r and so, due to orthogonality of the radial wave functions RnL{r), cannot connect the 
components originating from different (nL)-subshells. Instead, the quadrupole operator 
r^l2i does connect the components with different (nL). So, the scissors operators are more 
selective than the E21 operator. Every A'^ = 1^ state responds to both Ml and E21 external 
fields. Magnitudes of the responses depend on the wave function of the state. 
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The third and forth plots show STIRAP populations when the reaction runs through 
Xjx — 10 and Xjx — ll branches of the dipole plasmon (RF ElO and Ell cases, respectively). 
Like in the photoabsorptions, = 21 mode dominates over 20 and 22. In RF ElO case, 
the populations of both [211]-[202] and [211]-[200] leh states arc equally strong. Detection 
of these states provides energies of electron levels [202] and [200]. Besides, the scissors mode 
associated with the state [211]-202] can be observed. In RF Ell case, only [211]-[200] leh 
state is well populated. The differences between RF ElO and Ell cases can be of use (for 
example, to distinguish [211]-[202] and [211]-[200] leh states). 

B. Heavy oblate supported clusters 

The low-energy spectrum in heavy clusters is expected to be involved. In general, this is 
indeed the case. However, as is seen from Fig. 4 (second and third columns), the picture 
essentially depends on the reaction. Some responses reveal pronounced structures. We will 
analyse these structures by using Fig. 5 as a guide. Though Fig. 5 shows the level scheme 
for the small cluster, it hints the principle trends. 

The first line of Fig. 4 displays the quadrupole photoabsorption. We see bunching of the 
strength at 0.5-1.0 eV. The bunching involves two kinds of quadrupole transitions: i) between 
neighbouring subshells (like E20([220]-[200]) in Fig. 5) and ii) inside the valence subshell 
(hke E22([220]-[202]) in Fig. 5). The second line of Fig. 4 exposes photoabsorption for the 
scissors mode. The scissors peak rises and blue-shifts with the deformation. Unfortunately, 
both LEQM and scissors modes cannot be detected in photoabsorption and their plots are 
given here mainly for better understanding of the subsequent STIRAP results. 

The RF ElO plots demonstrate the overwhelming dominance of A// = 21 mode over 
20. This can be explained by statistical arguments. First, the mode A/x = 21 covers two 
projections = ±1 instead of one in Xji — 20. Second, the mode X^ — 21 usually has much 
more transition single-particle matrix elements than 20. 

The RF ElO response is presented by two distinct peaks. The first one originates from 
the deformation splitting (like E21([211]-[202]) in Fig. 5) and exhibits a blue-shift with 
increasing \52\- This peak is not seen in E2 photoabsorption. At the same time, it displays 
the correlation with the scissors mode. The second peak at ~ 0.6 — 0.7 eV is mainly 
determined by E21 transitions between neighbouring subshells (like E21([211]-[200]) in Fig. 
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5) and its energy does not depend on the deformation. 

In RF Ell case, the mode 22 comes to play. The response is entangled at 82 = —0.3 
but exhibits a pronounced low-energy (~ 0.1 eV) Xfi =22 structure at 62 = —0.5. For 
the first glance, this structure looks puzzling since it cannot be explained by any option 
discussed above. Indeed, its energy does not visibly depend on the deformation splitting 
(like E22([220]-[202] in Fig. 5 or 22 bump in Fig. 4, first line). It also has too low energy to 
be explained by transitions between the neighbouring subshells. Our analysis has revealed 
very specific origin of this structure. It is explained in Fig. 6 as a result of bunching [A^ri^A] 
levels with the same Uz (see E22 transitions [523] [521] and [532] — > [530]). The bunching 
signifies that in large systems with strong deformation the asymptotic Nilsson-Clemenger 
quantum number (number of quants along the symmetry axis z) becomes exact. 



C. STIRAP discussion 



Some essential points concerning detection of LEQM in STIRAP should be commented. 

The lifetime of the dipole plasmon is mainly determined by the Landau fragmentation 
while the radiative decay of the plasmon usually plays a minor role. The advantage of the 
STIRAP is that it strongly enhances the population of the desirable low-energy states in 
spite of a weak natural radiative decay to them from intermediate states. 

To distinguish in STIRAP the electron LEQM from possible excitations of other na- 
ture (contributions of impurities, ionic modes, etc), one can use the feature of electronic 
excitations to exhibit the inverse dependence of their energy on cluster size. 



It worth noting that STIRAP deals with both free and supported systems 



121. 



IV. CONCLUSIONS 



The properties of the low-energy quadrupole modes (LEQM) in deformed clusters were 
analysed within the Kohn-Sham LDA RPA approach. The stimulated Raman adiabatic pas- 
sage (STIRAP) was considered as the most suitable experimental method to detect LEQM. 
The population of the LEQM in STIRAP was calculated for two special cases: free light 
deformed clusters and supported heavy strongly-oblate clusters. In the first case, LEQM 
spectrum is dilute and can be resolved in STIRAP. The spectrum is easily associated with 
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particular electron-hole pairs. This finally allows to determine the single-particle energies 
near the Fermi level. Besides, close connection between the Xji —21 and scissors modes 
allows to observe the latter. 

In heavy clusters, STIRAP population through the AyU = 10 branch of the dipolc plasmon 
is strictly dominated by two A/i =21 peaks. The nature of the peaks is explained. In the 
population through the Xfj, — 11 branch, the soft X/j, —22 mode devotes a special attention. 
It originates from the fact that in large systems with strong deformation the asymptotic 
Nilsson-Clemenger quantum number becomes exact. 

It worth noting, that the LEQM properties described above are of a general character and 
so should take place not only in sodium clusters considered here but also in other clusters 
with quantum shells, for example in clusters of noble metals. Since any other information 
on electron LEQM in deformed clusters is absent, our study can be used as a first guide in 
the field. 
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Figure captions: 



Figure 1: Scheme of population of tlie low-energy quadrupole states =20, 21 and 
22 in RS, RF and STIRAP reactions through A/i = 10 (left) and 11 (right) branches of the 
dipole plasmon. 

Figure 2: Photoabsorption cross section for the dipole plasmon in Na^g, Na^g and Najy. 
Parameters of quadrupole and hexadecapole deformations are given in boxes. The exper- 
imental data 18[ (triangles) are compared with our results given as bars (for every RPA 
state) and the strength function smoothed by the Lorentz weight (with averaging parameter 
0.25 eV). Contributions of A/i = 10 and 11 dipole branches (the latter is twice stronger) are 
given by dashed curves. The bars are given in eVA^. 

Figure 3: The dipole plasmon in Nans at oblate deformations 82 =-0.1 (first line), 82 =- 
0.3 (second line), and 62 =-0.5 (third line). Contributions of A/i = 10 and 11 dipole branches 
(the latter is twice stronger) are given by dashed curves. The photoabsorption is smoothed 
by the Lorentz weight with A = 0.25 eV. 

Figure 4: Low-energy electron modes in Na^g (left column ) and Nang at 62 =-0.3 
(middle column) and 62 =-0.5 (right column). The plots exhibit LEQM photoabsorption 
(first line), scissors Ml photoabsorption (second line), STIRAP population of LEQM through 
A/i = 10 (third line) and 11 (forth line) dipole branches. LEQM are depicted by solid 
(A/i = 20), dashed (A/i = 21), and dotted (A// = 22) curves. All the responses are smoothed 
by the Lorentz weight with A = 0.1 eV. 

Figure 5: The electron level scheme for Na^g in the spherical limit (left) and at the 
equilibrium deformation (right). The Fermi level is [211]. Arrows depict the possible low- 
energy quadrupole hole-electron E2fi transitions. 

Figure 6: Deformation splitting of the electron subshell 2/ in Nang. The oblate defor- 
mations are listed below the plot. The levels are marked by Nilsson-Clemenger quantum 
numbers [NuzA]. The plot demonstrates the origin of the soft E22 transitions. 
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